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^vj Abstract. We present a general survey of heliocentric LISA orbits, hoping it 

, might help in the exercise of rescoping the mission. We try to semi-analytically 

1-^ optimize the orbital parameters in order to minimize the disturbances coming 

^_. from the Earth-LISA interaction. In a set of numerical simulations we include non- 

autonomous perturbations and provide an estimate of Doppler shift and breathing 
\I as a function of the trailing angle. 
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1. Introduction 



The LISA space experiment to detect low frequency gravitational waves has been for 
a long time a priority mission of space agencies, both in Europe and in the US. There 
has recently been an ample discussion on a possible scaled-down version of the LISA 
1^ mission that, in order to meet tighter budget constraints, could be characterized by 

00 a shorter arm length L, a closer mean distance from the Earth (a smaller trailing 

Cn angle) and maybe a 2-arins (4-links) configuration, giving up the third arm. In this 

[ — case, it becomes natural to consider a right angle geometry as an alternative to the 

^D traditional, 60°, equilateral triangle. 



Although other configurations are being evaluated by the ESA Concurrent Design 
Facility team, these triangular 'constellations' on heliocentric Earth-trailing orbits 
still remain the favorite choice. We focus our attention on the evaluation of the 
usual kinematic indicators of performance (arm flexing, breathing angles and Doppler 
shifts) when reducing both the size of the triangle and the Earth-LISA distance over 
the entire mission lifetime. As is well known, the interaction of LISA with the Earth 
is the major perturbation. The dominant effect is a parabolic drift characterized by 
a "rendez-vous" (RV) at which the distance between the constellation and the Earth 
is minimum. We investigate how additional perturbing effects influence the motion of 
LISA around the RV and how it is possible to optimize it. 

We assume the following guidelines: 

• Arm Length: L=l Gm. We consider two configurations: the equilateral triangle 
(ET) with side L and an isosceles right triangle (IRT) with two equal arms of 
length L and the third one L^/2 long. 
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• Flexing: the spacecraft (S/C) relative velocity in the sensitive axis (rate of 
change of the arm length) causes a Doppler shift of the laser frequency. We 
set a maximum bandwidth of 20 MHz over a 1/im carrier, corresponding to 
Vjc < 6.5 X 10~^, or T^ < 20 m/s. Note, for the sake of plot readability, that 
the Doppler shift (in MHz) has the same numerical value than the longitudinal 
velocity (in m/s). 

• Breathing angle: The relative motion of S/C's also imposes a continuous 
adjustment of the angle between two beams departing from the same corner, 
in order to track the opposite spacecrafts; this fluctuation over the nominal angle 
(60° or 90°) is referred to as breathing angle (BA). We demand BA < ±1.5°. 

• Trailing angle (TA, also referred to as Lag Angle): LISA follows the Earth on a 
heliocentric circular orbit and TA is the angle between the constellation and the 
Earth as seen from the Sun. TA is a good indicator of the Earth-LISA distance, 
because the radial secular motion (away from the Sun) of the constellation is 
normally much smaller than the tangential (along the Earth orbit) one. We 
demand TA as small as possible, compatibly with the above requirements. In any 
case, initial conditions are chosen in such a way that over the mission lifetime, 
TA never exceeds 21°. 

• Mission lifetime: 6 years. 

The plan of the paper is as follows: we start recalling the simple models describing the 
interaction between Sun and LISA (Sectionl2]) and Sun, Earth and LISA (Sectionlsl). In 
Section [4] we describe an optimization method with the aim of an important reduction 
of flexing, breathing angles and Doppler shifts. We test this method first on the 
simplified Sun-Earth-LISA model and then on a more complete model including the 
real gravitational effects due to the dynamics of the Solar System. Finally, in Section[5] 
conclusions are drawn. 

While graphs and details are given below, we anticipate here some results: 

• As far as Doppler and breathing requirements are concerned, a short LISA can 
safely be put in an orbit much closer to Earth: TA « 12° at RV (the baseline 
design was 20°), or 31 Gm. 

• Should we give up the third arm (keep only 4 optical links), a right angled triangle 
can be employed and performs at least as well as the usual equilateral triangle in 
several of the tests (Doppler, breathing, etc.) we carried on. 

2. Keplerian orbits 

We describe the interaction of LISA with the Earth in the framework of the Hill- 
Clohessy- Wiltshire (HCW) system [I] |2] . In this and next section, we will assume the 
Sun at rest in an inertial reference frame. The origin of the HCW frame rotates around 
the Sun on a circular reference orbit of radius Rq and with orthogonal axes oriented 
as follows: x is directed radially opposite the Sun, y is in the direction tangent to the 
motion and z is perpendicular to the ecliptic. The time evolution of these orbits can 
be described with adequate precision using the post-epicyclic approximation in the 
HCW frame [31ll|. 



2.1. Zero order approximation 



Under the effect of the Sun only, at zeroth order, the equations of motion for a S/C 
in the rotating frame are 



X — 2ujy — 3uj X 
y + 2ujx ^ 0, 
z + uP'z = 0, 



0, 



(1) 



where cj = ^JGMq/R^ is the mean motion and the most general solution is a 
combination of an ellipse in the xy plane and an oscillation in the z direction 
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The natural choice, in order to avoid drifts and offsets, is to set yo — —2ujxq and 
Xq = ujyo/2, so that the trajectory is reduced to a combination of simple oscillations 
along the three axes. Moreover, for a rigid, polygonal constellation, the distance of 
the S/C from the origin must be constant, say h, so that we obtain [51 
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Thc orbit of one of the S/C's in the HCW frame is a circular motion with constant 
angular velocity w and radius h around the origin in a plane inclined of ±60° with 
respect to the xy (ecliptic) plane. A second S/C, describing the same path (with a 
certain delay), will be at a constant distance i, from the first one. For n such S/Cs, 
on the vertices of a regular polygon, their distance h from the origin and the relative 
phase delay are 

2sm(7r/n) n 

Finally, we put xq ~ 1/2 to be consistent with the notations of ^^ and the zero-order 
orbits turn out to be 
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The ET constellation is obtained with n — 3 and k = 1,2,3, while the IRT one 
corresponds to n = 4 (and k = 1,2,3), namely 
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2.2. First order approximation 

A first optimization is possible by changing the tilt angle of the constellation, i.e. the 
inclination of the triangle with respect to the ecliptic. Introducing a parameter 6i [3] 
such that: 

±60° + Si 
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where h = i/^/3 for the ET configuration and h = £/V2 for the IRT. It can be shown 
[3] that, due to the small eccentricity and inclination of the orbits, the solution given 
by the above zero and first order terms differs from the exact Keplerian solution by 
less than 0.03% making the method of analytical series expansion a useful basis for 
an analytical model of the motions of LISA. 

Choosing 61 — 5/8, the first order (Keplerian) flexing is minimized in both ET 
and IRT configurations, giving, with arm lengths of order 1 Gm, an extra angle 
of respectively 7' and 4'. Table [l] reports some orbit indicators (flexing, breathing 
angles, and Doppler shifts) relative to both IRT and ET configurations for 61 = and 
di = 5/8. The indicators A"*" and A~ represent the difference between the maximum 
and minimum(respectively) value and the nominal value of a given parameter, over 
the 6 years of the mission. 



3. The Earth eflfect 

We now include the perturbation due to the Earth. 

In the analytic approach the Earth is assumed at rest in the Hill frame. As shown in 

[B] , in the rotating frame, the global dynamics in the coupled fields of Sun and Earth 

is characterized by secular terms producing, in the long run, a drift away from the 

Earth: this is linear in time in the radial direction and quadratic in the tangential 

direction. 

In order to minimize this quadratic y drift, an intuitive strategy is to choose initial 
conditions such that LISA is a little further out at start, approaches the Earth, reaches 
its minimum distance at mid mission and departs after that [7]. However, other 
strategies can be devised that provide better performance of the constellation. An 
appreciable reduction of the flexing due to the Earth tidal field is in any case possible, 
over a limited time span, by suitable tuning of all orbital parameters [51 [HI HO]- In 
our analytical approach, in order to keep things simple, we still use three identical 
orbits (apart for relative phase shifts, see ((sl)) for the 3 S/Cs of the constellation and. 



Table 1: Change of the relevant orbit indicators (arm length, breathing, Doppler 
modulation) for the three S/C's, for both IRT and ET configurations, in the standard 
^1 = and modified Si = 5/8 inclination. For each indicator, nominal value, average 
and deviations A+ and A^ (see Section 2.2), relative to the nominal value over a 
mission lifetime of 6 years are shown. 
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in addition to the above specs, we try and vary the tilt-angle Si and a subset of the 
initial conditions. We first consider a simplified model where the Earth describes a 
circular orbit of radius Rq — lAU around the Sun on the xy plane [HllS]. Introducing 
the minimum trailing angle TAq, taking place at time tg, the Earth coordinates 
(x®, y® , Zq) in the HCW frame are 

a;e = -Ra{l - cos(TAo)), y© = i?osin(TAo), z© ^ 0. 

We consider its effect as a constant + linear force to be added to the equations 
of motion. We define the distance of the Earth from the origin, d© = \/x'^ + y"^ , and 
introduce the perturbation parameter 

£© = ^1-K ^^..^ (6) 
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so that the effect of the Earth is given by the perturbation 

hx = e(^uj'^{x(^ + Ciix + Ci2y) 
hy = £(s^^{y(B + Ci2a; + C22y) 

fez = £(B^ (^e — z), 
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The equations of motion can be solved with the perturbation method and the solutions, 
to be added to the zero and first order solutions, are in the form 

xl '(t) = {Ak^x + Bk,xt) sinut + {Ck,x + Dk,xt) cosut + Ek^x + Fxt, 

yi^\t) = {Ak,y + Bk.yt) sinwt + [Ck.y + Dk^yt) cosujt + Ek,y + Fk,yt + Gyt^, (8) 

^k (i) = i^k,z + Bk,zt) sin ut + {Ck,z + Dk,zt) cos ujt. 
The integration constants Ai ■ ■ ■ Gy are defined by the choice of initial conditions (see 



Appendix A). The secular terms appearing in the solution generate the parabolic drift 



around the RV (see Appendix B I . 

By collecting terms, the orbit of the S/Ck is 

r,(f)=ri°)(t)+r«(t) + rf)(i) (9) 

where the zero and first order terms are respectively given by (|3|) and (Wl) and ri, (t) 
is given by (|8|). The terms growing as t and t^ in (respectively) x^-^^(t) and y^^\t) 



vanish when one calculates the relative motion between S/Cs, being Fx and Gy equal 
for all S/Cs, therefore the increase in fiexing with time is only due to secular terms as 
isini and tcost. 

It is now useful to define some quantities in the heliocentric frame: 

• unit vectors of the rotating frame axes: 
Ux = {cos(wt - TAo), sin(wi - TAq), 0}, 
Uy = {-sin(a;i- TAq), cos(wt - TAq), 0}, 
u. = {0,0,1} 

• Position of the Earth: R©(t) = Ro{cosujt,smujt,0}. 

• Orbit of S/Cfe: 

R/c(0 = (i?o + xk{t))ux + yk(t)uy + zk{t)uz (10) 

• LISA barycenter position: Rg(i) = | J2k '^k{t)- 

• LISA arm vectors: Rij(t) = ^ij{t) ~ ^j{t) ^ ri(^)- 

• LISA arm lengths: i,j(t) = |ry(f)| = |Rj(t) - Rj(t)|. 

• Doppler shifts: Vij{t) = j^Lij{t) 

• Trailing angle: 



TA(0 = l^arccosf^#-^ 
• Breathing angles: 

e,[t)- ^ arccos(^^^^_^^^^^_^^^^ 



4. Minimization of the flexing 

We shall address the choice of an orbit that minimizes the flexing of the arms in three 
steps: first, by optimizing with respect to the tilt angle only. Then, by perturbing 
the initial conditions of the three S/Cs in a still analytic approach. Finally, with a 
fully numerical integration of the S/C orbits, taking into account all major perturbing 
effects. The mission begins at the time tini and is assumed to last At ~ 6 years. We 
refer to "mid-mission" or tmid = tini + 3 yrs the time half way into the mission. 

In the panels of Figure ll] we consider both the ET (red) and the IRT (blue) 
configurations: for these plots the optimization is only done by evaluating the optimal 
tilt-angle over the mission lifetime of 6 years. Optimization is performed by minimizing 
the RMS flexing of the 3 arms over the entire mission duration. The time of closest 
approach to the Earth is 3 years after the mission starts and this identifies TAq. 




Figure 1: Effect of the Earth (non optimized): left: maximum Doppler shift in ni/s. 
Right: maximum change of the angle between L12 and L13 (short arms in the IRT). 
Red curves: ET, blue curves: IRT. 



The requirement of 6.5 x 10^^ on the Doppler shift shows that we can reduce the 
minimum TA to less than 10° for both cases (8° for the IRT). By choosing an optimal 
tilt angle we can strongly reduce the rms flexing of the arms; this angle always turns 
out to slightly differ from the canonical 60°. On the other hand, the breathing angle 
requirement of ±1.5° sets a limit at TA ~ 14°. Breathing appears therefore the major 
obstacle to an appreciable reduction of TA. In the next subsection, we attempt a more 
general optimization strategy. 



4.I. Cost function 

In order to extend the optimization to a wider set of parameters, we need to introduce 
a "cost function" , i.e. a suitable function of the relevant quantities (the flexings) that 
has to be minimized. To this purpose, we define the following cost function, suitable 
for both configurations (and therefore different from those proposed in [5]): 



= ((il2 - (^12))' + {L23 - {L2z)f + {L3I - {L3l)f) 



(11) 
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indicates average over the mission time. Although in the IRT case the third arm is 
not monitored and its flexing could appear as a useless burden to the cost function, 



we maintained the same cr^, as defined in ( |11[ ), for both configurations: minimizing 
the flexing of all arms is a way to render the triangle " more rigid" and is an effective 
strategy, as we shall show, to minimize the angle breathing as well. 

4-2. Perturbation of initial conditions - semi-analytic approach 

As we extend our optimization strategy, remaining as close as possible to an analytic 
approach, we must restrict the space of free parameters. Our choice is a subset of the 
initial conditions of the unperturbed orbits. The Earth produces linear and secular 
terms in the orbits, as shown by (l8|. On the other hand, (pi) show that, in general, 
a linear drift exists in the y component (that we canceled by setting 2loxo + j/q = 0) . 
We can therefore choose a suitable offset in the initial conditions in such a way that 
these two linear drifts compensate each other. 

Therefore, if we set 

xo.k^-^+ek, fc = l,2,3 (12) 

the first two components of position vector ([3]) become 

(0) ^ 1 

xl = - cos cTfc + efe(4-3cosa;i), 

2sm 7r/n) 2 , , 

y (13) 

Vk = r. ■ / ; N sin CTfc - 6efc(a;t - sinwt) 
2sm(7r/n) 

while the third one is unchanged. 

In this way trajectories and arm-lengths are affected by a perturbation which grows 
linearly in time. In the expansion (l9| we modify only the zero-order terms: in 
principle, the variation should be propagated through the higher-order terms, but 
the contribution is of order e£/R in the first-order terms and even smaller in the term 
describing the Earth effect: we can therefore safely neglect them. 

The required amount of this variation can be determined by minimizing the cost 



function ct (£1,62,63, (5i) defined in (11). However, the analytic expression for ct is 
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sufficiently cumbersome to impose a numerical minimization: this, on the other hand, 
allows us to use the exact equations of motion: 

Jx.k 

(14) 



where ffc is the Sun-|-Earth force per unit mass acting on the fc-th S/C, expressed in 
the HCW coordinate system. 

Setting the RV at tg = tmid, for the IRT and ET configurations respectively, the 
minima correspond to 

IRT : (5i = 0.808, ei = 867 km, 62 = 519 km, 63 = 66 km, , , 

15) 
ET : (5i= 0.894, ei = 523 km, 62 = 64 km, 63 = 7 km. 



Table 2: Variation of the same orbital indicators as in Table IT] (arm-length, 
breathing, Doppler modulation) including the Earth effect (assumed on a circular 



orbit) corresponding to the optimal data of ( 15 ) for the IRT and ET constellations 
(left and right, respectively). 
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The results of the optimization are reported in Table [2] The trailing angles in both 
cases at Uni and t/i„ are 12.8° degrees (33 Gm from the Earth). The improvement in 
the values of the performance indicators in the optimized cases is quite evident. 



4-3. Numerical optimization 

In this section we describe the fully numeric evaluation and minimization of the 
cost function ( |11[ ) by solving the exact equations of motion and taking into account 
perturbing effect of the Sun, Venus, Earth, Moon, Mars and Jupiter. Their real 
trajectories TZq (t) , T?.® {t),TZ^ (i) , TZc^ (t) , TZ^ (t) , TZq^ (t) in the Solar System Barycenter 
(SSB), are provided by the JPL HORIZON ephemerides [Hj, with the following 
characteristics: 

• Reference epoch: J2000.0 

• XY-p\a,ne: plane of the Earth's orbit at the reference epoch. 

• X-axis: out along ascending node of instantaneous plane of the Earth's orbit and 
the Earth's mean equator at the reference epoch. 
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• Z-axis: perpendicular to the XY-plane in the directional (+ or -) sense of Earth's 
north pole at the reference epoch. 

• step: 1 day. 

In the simplified model used till here, the Sun is assumed at rest at the center of an 
inertial frame. However, the true inertial frame is represented by the Solar System 
Barycenter (SSB), where the Sun moves in a non negligible and complex (non simply 
periodic) way: in this frame the motion of the Earth is substantially different from 
an ellipse, and therefore the initial condition that we adopted for the S/Cs using ^ 
are no longer suitable. Moreover, the motion of the Sun is also a relevant source of 
perturbation. Therefore, to account for these additional effects while maintaining the 
convenient, Sun-centered HCW description, we must complete the equations of motion 
with an apparent force deriving from the acceleration of the Sun relative to the SBB. 



The equations of motion are as (14), with the forcing term modified as follows 



/aj.fc = (ffc -7^0) -u^i;, /j^,fe = (ffc - 7^o) • u,y /^^fc = (ffc -7^0) -u 

where f^ is the total Newtonian force per unit mass on the fc-th S/C. 
GA/„ 



GM 



and Rfe is the position of fc-th S/C in the heliocentric frame (given by (lOl) 



The amplitude of flexing and breathing scales inversely with the LISA-Earth distance. 
This can be intuitively explained as follows: a small flexing is obtained if the 
constellation rapidly moves away from the Earth, its main source of perturbation 
to a rigid configuration. However, the overall distance in the mission lifetime must be 
bound within reasonable values dictated by communication requirements. 

An analytical study of the evolution of the Earth-LISA distance is shown in | Appendix] 
p3| were it is verified that the LISA-Earth distance increases as t^, after (and before) 
the RV. Moreover, there is an additional sinusoidal modulation at 1 year period due to 
the eccentricity of the Earth's orbit. We prove that the minima of the sinusoid occur 
at well defined epochs that depend on the allowed minimum TA but not on the epoch 
to of the RV. Therefore, in order to minimize the Earth-LISA distance, the optimal 
choice for ip is just one of these minima (B.5 1. This shows, as mentioned in Section^ 



that other choices of tg, different from tmid, can minimize flexing and breathing. 

In the following we shall discuss two cases: to = tini and io = imid- The value of 
TAq is chosen as the minimum one that allows a breathing angle smaller that 1.5°, as 
required. For the two conflgurations and the two kinds of RV considered the minima 
of the cost function is found at the following values of parameters; 

IRT - RV at the beginning, (to = tini): 

Si — 0.061, £1 — 430 fcm, e2 = — llSfcm, £3 — —9 km. 

IRT - RV at mid mission (to = tmirf): 

,5i = -0.290, ei = 28 km, £3 = -55 km, €3 == -170 km. 
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ET - RV at the beginning (to = ^ini): 

5i = 0.473, £1 = 70 km, e2 = -483 km, £3 = 42 km. 

ET - RV at mid mission (tg = tmid)'- 

5i = -0.047 £1 = 193 km, £2 = 52 fcm, £3 = 18 km. 

Table Is] provides more results and details for the four cases (2 configurations x 2 RV 
times) considered here. 

Some results are also plotted in Figure [2] and |3] for IRT and ET configurations, 
respectively. The ranges of LISA-Earth distances and trailing angles, as well as the 
initial conditions for the S/Cs are reported in Table l4] We observe that the minimum 
TA is larger when the RV is at mid mission, but the ATA is smaller. In general, 
breathing angles within the specs of 1.5°, can be obtained at smaller distance from 
Earth for the ET than the IRT configuration. 

Figure l4] shows, in a synoptic way, the results of our optimization procedure with 
respect to Doppler and breathing angle, vs. the minimum trailing angle TAq. By 
comparing these optimized results with those derived from the simplest model shown in 
Figurefll we see that, even considering many more perturbing agents, the optimization 
manages to reduce both performance indicators by about a factor of 2 at small TA. 
Again we see that the requirement on the breathing remains the most stringent 
constraint. However, while complying with keeping the breathing within ±1.5°, we 
can address the reduction of flexing following again two opposite strategies: we can set 
the RV at the beginning of the mission, achieving the lowest values of TAq (we have 
12.5° for the IRT and 12.1° for the ET), and accept a maximum ATA of about 8.5 
degrees in both cases. Else, if RV takes place at mid-mission, we must accept larger 
values of TAq (13.9° for the IRT and 13.8° for the ET), but TA will change much less 
during the mission: ATA is less than 4.8 degrees in both cases. 

5. Conclusions 

We have shown that the choice of heliocentric orbits for LISA is a viable solution even 
when reducing the arm-length: this allows a substantial reduction in the TA (with 
deriving beneficial savings for placement in orbit and communications with Earth) , of 
an amount that depends on the assumed mission duration. For an expected mission 
time of 6 years, the minimum value of the TA can be reduced to about 12°. Should a 
2-link interferometer be preferred for a new, cheaper version of the LISA mission, the 
Isosceles Right Triangle is a viable configuration, as stable as the Equilateral Triangle 
in all of the tests we have computed. The amount of flexing that the constellation 
undergoes during the mission depends strongly on the initial conditions. The reasons of 
this behaviour lie mostly in the time dependent perturbations due to the eccentricity 
of the Earth orbit and to the motion of the Sun with respect to the Solar System 
Barycenter. A more detailed analysis of these effects is underway. 
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Figure 2: Fully numerical optimization (Section 4.3 1 with respect to both initial 
conditions and tilt angle for the IRT configuration. Left panels: RV at the beginning 
of the mission, right panels: RV at mid-mission. Top panels: breathing angles. Center 
panels: Doppler shifts. Lower panels: distance LISA-Earth, expressed as TA (red lines 



are obtained using (B.6) 
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Figure 3: Same as in Figure |2J but for the ET configuration. Left panels: RV at the 
beginning of the mission, right panels: RV at mid-mission. Top panels: breathing 
angles. Center panels: Doppler shifts. Lower panels: distance LISA-Earth, expressed 



as TA (red lines are obtained using (B.6) 
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Figure 4: Results of the optimization. Left: maximum Doppler shift in mj s. Right: 
maximum change of the angle between L\2 and I/13 (short arms in the IRT). Red 
curves: ET (RV at the beginning of the mission), blue curves: IRT (RV at the 
beginning of the mission), yellow curves: ET (RV at mid mission), green curves: 
IRT (RV at mid mission). 

Table 3: Variation of the same orbital indicators as in Table 111 and p] (arm length, 
breathing, Doppler modulation) including the effect of the main bodies of the Solar 
System for the IRT and ET constellations (left and right, respectively). 
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Table 4: Initial conditions for the IRT (top) and ET (bottom) configuration in the 
heliocentric reference frame. 





Xito) 
[Cm] 


Yito) 

[km] 


Z{to) Xito) 
[km] [km/h] 


Y{to) 
[km/h] 


Zito) 
[km/h] 


IRT, RV at tini 
S/Cl 

S/C2 
S/C3 


149932288 
149588968 
149222635 


t,nin= 2018-Oct-05 
2412779 612461 1721 
1697856 2895 1471 
2401359 609567 1729 


106957 
107213 
107465 




438 




IRT, RV at t„nd 
S/Cl 

S/C2 
S/C3 


149753118 
149404775 
149042293 


trnin= 2018-Oct-07 

2563623 611220 1831 
1872191 11110 1595 
2568312 608484 1845 


107004 
107265 
107514 


10 

438 

9 


ET, RV at tini 
S/Cl 
S/C2 
S/C3 


149884804 
149453230 
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t„,in= 2018-Oct-05 
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51258 248057 217 
1052373 248057 576 
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311 


ET, RV at trnid 

S/Cl 
S/C2 
S/C3 


149701044 
149267830 
149266254 


t„,i„=2018-Oct-07 
1382500 499337 987 
902343 257543 824 
1881003 237763 1167 


107064 
107378 
107370 


9 
307 
314 
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Appendix A. Initial Conditions for the motion in the Earth field. 

Equation (l8| can be recast in the following, equivalent but more explicit form, where 
the constants Cy of ^ are folded into the solution: 

(2) 

X f, = X(^ + 2Ak + Sy^wt + Bk cos ujt + Ck sin ujt + 

(Cii + I2C22) cos cTfc + 2 (2Ci2 + (Cii +4C22)c^t)singfc 

yf = 42/e - (3^ + 2xe)wt - -ye(wi)' + Dk 
+ 2{Ck cos tot — Bk sin ujt) 

(aCii + I6C22) cos gfc ~ 2(2Ci2 + (Cii + 4C22)c^^) cos gfc 

4^/3 

(2) ^ 

zl' ^ Ek COS CTfc + Ffc sm CTfc - jCFk sm CTfc , 

where the CTfc are the time-dependent phases ([3]). 

The 18 constants Ak, ...,-Ffc are determined by the initial conditions, that are chosen 
assuming (xf \ y[^\ 4^^)=(0,0,0) at i = 0. 

They are: 

. _ „ 2C22 cos CTfc - Ci 2 sin CTfc 

^'-^^ 71 ' 

Bk = - Xfs - ^^{{Cii-4C22) cos <tI+4Ci2 Sinai), 

./•if 
Ck = -2ye + ^((Cn-4C22)sinCT^-4Ci2COSCT^), 

Dk = ^[Ci2Cosal-2[Ci^+iC22) Sinai), 
Ek^ -^ sin^ CT°, ^" = -1 (^'^^ " '^"2ct°) , 

where the CTfc — —^^ — - are the relative phase shifts of Jsfc evaluated at t = 0. 

Appendix B. Distance Earth-LISA barycenter in epicyclic approximation 

Here we calculate the analytic expression for the distance between a particle (i.e.: the 
LISA barycenter) and the Earth, taking into account the eccentricity of the orbit. 

We consider the Earth orbit in epicyclic approximation (at t = in the perihelion) in 
the inertial frame centered in the Sun: 

R©(i) = Rolcosut + - cos 2ujt — -e, sinwt + - sin2a;i, 0}. (B.l) 

The LISA barycenter, as a first approximation, can be considered at rest in the HCW 
frame at TAq degrees from the Earth. In the inertial frame its trajectory is 

-Rg^t) = i?o{cos(a;t - TAo),sin(a;t - TAo),0} 
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At zero-order, the force of the Earth on the particle is 

f = eRoui^lfx + ecx cos cut + es^ sinuit, fy + ecy cos cut + esy sinojt, 0} (B.2) 

where e ss 0.01671 is the eccentricity of the Earth's orbit and 

Afffi 2 

U = -2: fy = 



4M0V2-2cosTAo' ' ^ tanTAo/2' 



2 1 

1: c„ 



cosTAo-1 "' ''^ tanTAo/2' 
« -2. 



tanTAo/2' " cosTAo-1 

The new perturbation parameter (e = 2.6 x 10"** for TAq = 10°) is shghtly different 
from that introduced in ^ to show the exphcit dependence on TAq. The coefficients 
eCx-, eCy, esxi eSy are much smaUer than fx and fy and therefore we neglect the terms 
proportional to ee in ( |B.2[ ) and solve perturbatively the HCW equations, assuming 
Vg = {0,0,0} as the unperturbed motion. We calculate the perturbation ri(i) with 
the assumptions ri(io) = {0, 0, 0} and ri(io) = {0, 0, 0} where to is the epoch at which 
we put the particle at TAq degrees from the Earth. Letting t' = t — to we have 

ri(i') = eRoifxil - coswi') + 2fy{ujt' - smcut'), 

2fxisiiiu}t' - ujt') + fy{4: - 4:C0SUjt' - 3/2a;2t'2), (B.3) 

0} 



We transform (B.3 1 in the inertial coordinates using ( 10 ) and we calculate the distance 
d(t) of the particle from the Earth using the expression (B.l). Finally, we expand in 
Taylor series the distance to the first order in e and e 

d{t) = do + edi{t) + ed2{t) + 0(e2) (B.4) 

where 



do = i?oV2-2cosTAo 

di{t) = ° === [(cosTAo - l)cosa;i + 2sinTAosinwt] 

V 2 — 2 cos TAq 

dM')= , ^° y 

^ ^ 2V2-2cosTAo 

X fx + 2fyijt' - fx cos TAo - 2fyOJt' cos TAo + l/2(-8^ + 4/,tjt'+ 
+3fyLu'^t''^) sin TAo + cos ujt'{-fx + fx cos TAo + 4^ sin TAo) + 
+ i~2fy + 2fy cos TAo - 2fx sin TAo) sinujt' . 

The term do is a constant, the term di is a sum of sinusoids with lyear period. The 
(^2 term contains linear and quadratic terms in t — to, and is therefore negligible for 
t « to because e <^ e but it becomes dominant for larger t. 

The epochs of the relative minima and maxima of d{t) depend on TAo but not on to. 
They are found by equating to zero the first derivative of di : 

2 cos ojt sin TAq + (1 — cos TAo) sin uit = 0. 
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With the additional condition on the second derivative 

(1 — cos TAq) cos iOtmin — 2 sin TAq sin UJtmin > 0, 

the minima occur at 



1 



rain^k 



arctan 



tan(TAo/2) 



2fc7r 

1 



ke 



In the same fashion of (B.4), the TA can be obtained, to first order in e, as 



TA(t) = TAo + 2e sin wi + e 



4/y(cosa;t' - 1) + 2f^{ujt' - smcot') + -fyitot'f 



(B.5) 



.(B. 



Although the epochs tmin are obtained using a first-order approximation, they are 
in good agreement with the exact values (for an example, see the bottom panels of 
Figure [2] and |3|. 
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